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Transverse spin diffusion in a polarized, interacting Fermi gas leads to the Leggett-Rice effect, 
where the spin current processes around the local magnetization. With a spin-echo sequence both 
the transverse diffusivity and the spin-rotation parameter 7 are obtained; the sign of 7 reveals 
the repulsive or attractive character of the effective interaction. In a trapped Fermi gas the spin 
diffusion equations become nonlinear, and their numerical solution exhibits an inhomogeneous spin 
state even at the spin echo time. While the microscopic diffusivity and 7 increase at weak coupling, 
their apparent values inferred from the trap-averaged magnetization saturate in agreement with a 
recent experiment for a dilute ultracold Fermi gas. 

PACS numbers: 67.85.Lm, 75.76.+j, 51.10.+y 


I. INTRODUCTION 

Transverse spin diffusion occurs when the magnetiza¬ 
tion is oriented along different directions, for instance, in 
a spin helix. It has been observed in spin-echo experi¬ 
ments in helium [I|, polarized hydrogen, and, recently, 
ultracold atomic gases M- The transverse magnetiza¬ 
tion evolves according to a diffusion equation, but there 
are additional terms from the precession of the spin cur¬ 
rents around the local magnetization. This Leggett-Rice 
effect @ is related to the identical spin-rotation effect 
Q and leads to reactive spin currents and spin waves, 
which have been observed in ultracold Fermi gases Q- 
Consider a polarized sample with small transverse mag¬ 
netization (small tipping angle). The transverse magne¬ 
tization is conveniently combined into a complex number 
TO+ = uix + ivTiy , which evolves with a complex diffusion 
coefficient, 

dm+ ^2 • 

^-V m+ - iaa;3m+, (1) 

ot 1 — i^ruz 

where is the transverse diffusivity and /.tm^ = 7 (at 
full polarization) the dimensionless spin-rotation param¬ 
eter. For small p. —?> 0 this is the usual diffusion equation, 
while for large |^| the diffusion equation has an imaginary 
effective diffusivity and resembles the Schrddinger equa¬ 
tion In a spin-echo pulse sequence, the second term 
in Eq. m expresses a linear gradient a of the external 
magnetic field along the X3 direction, which winds the 
magnetization into a helix. After a time a tt pulse is 
applied, which is equivalent to flipping the sign of a. In 
the ensuing time evolution, the helix unwinds until the 
magnetization is realigned at the echo time te = 2t^. In 
the presence of spin rotation /r yf: 0, realignment at 
occurs at a phase angle (j) rx yi with respect to the ini¬ 
tial orientation at t = 0 and the value of /r can be 
inferred. For a strongly interacting Fermi gas the value 
of fj, has recently been measured and used to determine 
the spin-antisymmetric Fermi-liquid parameter F^, while 
the sign of fj, reveals the attractive or repulsive character 
of the effective interaction [^. Theoretically, Dq and /i 


for dilute, homogeneous Fermi gases have been computed 
using kinetic theory iH- 

In an infinite homogeneous system where the only po¬ 
sition dependence arises from the magnetic field gradient 
a, the phase angle (j) is directly proportional to the mi¬ 
croscopic parameter /i. Instead, for a finite homogeneous 
box the phase (j), and hence the apparent value of /r, sat¬ 
urates when the system size is reached [HI- Experiments 
with ultracold atomic gases typically employ a harmonic 
trapping potential which is both finite and inhomoge¬ 
neous: in this case, the diffusivity Dq, the Leggett-Rice 
parameter yt, and the magnetization ruz are strongly po¬ 
sition dependent, and the diffusion equation m becomes 
nonlinear. In previous studies the evolution has been 
linearized in order to determine the collective mode fre¬ 
quencies and decay rates in the trapping potential [l^ . 
Other studies consider the nonlinear evolution of the full 
phase-space distribution for a nondegenerate or collision¬ 
less trapped gas [l^. Here, I numerically solve the non¬ 
linear evolution equation with position dependent kinetic 
coefficients including medium scattering to obtain the 
transverse magnetization decay and the growth of the 
phase (j) for the specihc trap geometry used in the exper¬ 
iment. 

Kinetic theory 0 is employed to compute the spin 
evolution of the trapped gas. This method is well con¬ 
trolled in the weak-coupling limit, which is also the pa¬ 
rameter regime where finite-size corrections due to the 
trapping potential are most pronounced Q . On the other 
hand, at strong coupling toward unitarity and at low 
temperature near the superfluid phase transition of an 
attractive Fermi gas [T^ , kinetic theory is expected to re¬ 
ceive quantitative corrections from the effects of pairing 
and short quasiparticle lifetimes, which are incorporated 
for instance when computing transport from the Kubo 
formula within a Luttinger-Ward approach [110. 

This paper is structured as follows: in Sec. El the spin 
evolution equations for a trapped Fermi gas are derived 
from kinetic theory, while Sec. IIIII compares them to the 
known homogeneous limit. Section IIVI presents the re¬ 
sults for the magnetization profiles and apparent diffu- 
sivities, and Sec. IVl concludes with a discussion. 
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II. KINETIC THEORY 


Transport in an interacting Fermi gas may be described 
by kinetic theory for quasiparticles, as long as they are 
sufficiently long-lived. For a multi-component Fermi gas 
with two or more spin species one has to compute the 
time evolution of the spin distribution which is a 
matrix with components = npaa'{x^t) in spin space. 
The evolution equation derived by Landau and Silin [l^, 
[T^ reads 


dn 1 1 

+ 2 ^rnp]+ - 2 [VrSp 


-i 


+ -[£p,np\- = 


drip 


-W 

r coll 


where = Spaa'ix^t) is the matrix of single-particle 
energies. The left-hand side constitutes the drift term, 
while the right-hand side describes the change in the dis¬ 
tribution caused by collisions. Specifically for the spin- 
1/2 case. Up and Sp are 2x2 matrices in spin space, for 
instance, in the f, i basis. 

The spin matrices can be decomposed in terms of the 
identity / and the Pauli matrices ct. The occupation 
number matrix is written as 


Rp — 2 ’ 


■R-) 


( 3 ) 


where fp{x,t) is the particle number distribution func¬ 
tion and crp{x,t) the spin vector distribution in Bloch 
space. Similarly, the energy matrix 


Sp = epl + hp ■ a ( 4 ) 

combines the spin-independent single-particle energies 
ep{x,t) and a magnetic field hp{x,t). One may then 
rewrite Eq. m as 


_ dsp dfp 
dpj dxj dxj dpj 

^ dhp dcTp dhp dcTp 

dpj dxj dxj dpj 

and 




I 

dt 


E 


dsp dcTp 

dpj dxj 


dsp dcTp 

dxj dpj 


dl^d^ _ dhp dfp ] _ 2 

dpj dxj dxj dpj_ h ^ ‘ 


( 

V dt /coll 


( 6 ) 


The spin-rotation term hp x cTp is responsible for the 
Leggett-Rice effect. The single-particle energies are 


for a Fermi gas in a harmonic trapping potential V{x), 
which can be anisotropic with different trapping frequen¬ 
cies ujj in spatial direction j. The magnetic held 

hp{x,t) =-^fi{x,t), n = no + fimf, ( 8 ) 

w 

floix,t) = a{t)x3Z, Q.^i=—m{x,t) ( 9 ) 

is written in terms of the Larmor frequency = Oq -f 
rimf- For a spin-echo protocol it has two contributions: 
(i) flo = ^B{x,t) is due to the external magnetic held 
B {x,t) with gyromagnetic ratio 7. A spatially constant 
B is compensated by going to the co-rotating frame in 
Bloch space, but a magnetic held gradient of slope a 
along the X3 direction winds up the local magnetization 
into a spin spiral, (ii) The second contribution to the 
Larmor frequency, rimf, is a mean-held term proportional 
to the local magnetization m of a polarized Fermi gas. 
It leads to the precession of the spin current around m. 

The evolution of the full distribution functions fp{x, t) 
and crp(x,t) is simplihed by considering moments with 


respect to momentum: 

m{x,t) = j R^eTp{x,t) (12) 


with bare velocity Vpj = dsp/dpj = pj/m. An additional 
contribution fpdhp/dpj to the spin current is absent for 
momentum-independent hp. The local polarization is de- 
hned as M{x,t) = m{x,t)/n{x,t) where \M{x,t)\ < 1 . 
The spin current Jj is both a vector in Bloch space (bold 
symbol) and a vector in position space (j index): it en¬ 
codes how the magnetization changes as one goes along 
the j direction. The evolution equations for the moments 
read, using the specihc form of the single-particle energies 

atn + ^V,J/ = 0 (14) 

J” -1- anVjU + uj'^XjU =-A (15) 

'^n 

dtm + ^ Vj Jj -I- rio X m = 0 (16) 

j 

dtJj -I- a||P|| VjTn -I- a±P±\7jm + ojjXjfn 

+ ('’» +T”*) = (^i.„ ■ ™ 


£p{x,t) = ^ + V{x), V{x) = ( 7 ) 

j 


The projectors Pya = (a • m)rh and Pj_ = 1 — Py give 
the component of the magnetization gradient parallel and 
perpendicular to the local magnetization, respectively. 
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When deriving these equations, the mean field ilmf has 
been retained only in the spin-rotation term and not in 
the Poisson brackets in Eqs. m and This can be 
justified as the leading order in a controlled large-iV ex¬ 
pansion [ 2 ^. In general, the evolution of the currents 
depends on the second moments of fp and (Tp which 
in turn depend on higher moments. However, near lo¬ 
cal equilibrium the linearized Boltzmann equation relates 
the higher moments to the lower ones via the coefficients 
an, Tn, <a||,_L 5 td) and W, and one obtains a closed set of 
evolution equations. These coefficients are discussed in 
Sec, nm 


A. Initial conditions 


In the absence of the external magnetic field gradient 
a, the local equilibrium Fermi distribution is 

np±{x) = - ] -—— (18) 

exp{(3{ep{x) - /r±)) -h 1 

in terms of the chemical potential fi± of the majority 
(minority) component. The resulting density profile of 
the fully polarized gas is 

n{x) = Li 3 / 2 (-z+e“'^^^"’)) (19) 

with thermal wavelength A = {2'Kh '^and polylog¬ 
arithm Lis( 2 ;). The fugacity z+ = exp(,9/i+) of the major¬ 
ity component is determined by the total particle number 
N in the trap of average frequency w = (wia; 2 W 3 )^/^, 


N = - 


(/3&u)' 


■Li3(-2;+). 


( 20 ) 


At high temperatures in the Boltzmann regime the den¬ 
sity profile is Gaussian, n{x) oc exp(—/3B(a:)), while at 
low temperatures it approaches the Thomas-Fermi pro¬ 
file, 


n{x) = no(|l - ^ 
j 



3/2 


( 21 ) 


with Thomas-Fermi radius Rtfj = {2Ef and 
the Fermi energy Ep = {6Ny^^huj of a fully polarized 
gas. The density profile is shown in Fig.[T]both for weak 
coupling on the BCS side (upper panel) and for strong 
coupling at unitarity (lower panel). Note that the phase- 
space density X^n (solid black line) is unaffected by s- 
wave scattering in the fully polarized gas. 

In the spin-echo protocol Q, the gas is initially fully 
polarized, |m| = n or \M\ = 1 , at a tipping angle 6 from 
the z axis in Bloch space. 


m{x,t = 0) = {sm9,0,cos9)n{x). (22) 


This distribution is stationary for a = 0 and the particle 
and spin currents J" and Jj vanish. When the exter¬ 
nal gradient a is switched on along the x-i direction, the 




FIG. 1. (Color online) Kinetic coefficients and density pro¬ 
files in the trap on the BCS side (upper panel) and at unitar¬ 
ity (lower panel) at degenerate temperature (TlTp)i = 0 . 2 : 
total phase-space density X^n{x) (solid black line), diffu¬ 
sion parameter I3mai\ [dotted (brown) line], diffusive scat¬ 
tering rate tfItd [solid (red) line], interaction parameter 
tfX~^\W\ [dash-dotted (blue) line], and Leggett-Rice param¬ 
eter 7 = —ToWn/h [dashed (magenta) line]. 


distribution remains independent of xi and x^, but cur¬ 
rents are generated which can change the distribution 
along the x^ direction. Hence, the spin evolution effec¬ 
tively reduces to a one-dimensional problem for m and 
J 3 along the gradient direction. 


B. Kinetic coefficients 


The coefficients ay _L(a;) in Eq. (HB parametrize the 
strength of the spin current generated by a magnetization 
gradient. Their values are determined by the Boltzmann 
Muation linearized around the local equilibrium solution 
iP. There are different contributions from the longi¬ 
tudinal magnetization gradient due to the trap potential, 
and the transverse magnetization gradient from the helix. 
For a fully polarized gas. 


a\\{x) 

a±{x) 


an{x) = 

Pjx) 

mn{x) 


n{x) 

mx{x) 


(23) 

(24) 


where x(a:) = —A j3lj\i/2{—z+e is the local sus¬ 

ceptibility, and P{x) = —A“^/3“^ Li 5/2 (—the 
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local pressure of a free Fermi gas. In Fig. [1] is plot¬ 
ted as the dotted line: for a high temperature or low 
density in the outer regions of the trap it reaches the 
Boltzmann limit ay = a_L = 1/Pm, while for a high den¬ 
sity at the trap center it is enhanced and approaches the 
low-temperature limit an = Vp/3 la. In general, ay and 
a_L acquire interaction corrections [9j , but for large initial 
polarization and s-wave scattering these are negligible. 

Second, the coefficients Ty(a;) and T_L(a;) parametrize 
the decay of the spin current due to scattering Q , 


( 


dt 


coll 


d^p dsp / d(Tp\ 
(27r?i)3 dpj V dt 

^11 


T± 


(25) 


The component of the current P\\Jj parallel to the lo¬ 
cal magnetization m decays with the longitudinal diffu¬ 
sion time Ty, while the transverse component P±Jj = 
(1 — F)|)Jj decays with the transverse diffusion time tj_. 
Both differ for a polarized, strongly degenerate Fermi gas 
i 113, mi: all the states between the majority and the 
minority Fermi surfaces of a polarized gas are available 
for transverse scattering and tj_ can be much lower than 
Ty. For the experiment at hand Q, however, the low¬ 
est temperatures reached at the trap center are around 
T/Tp > 0.3, and previous studies have shown that in 
this temperature range ry and t± are nearly equal 0 - 
It is therefore jiKtified to work with a single decay time 

Td =Tj_Ki Ty 0 , 


h 

T_L 


sinh(/3/i) 1 


d^pi ... d^p4 


c± {2-Khf J 

X (5(3) (pi + P2 - P3 - Pi) + £p2 - 
X |T(pi,P 2)P [e“^''ni+n2+ -f e^'‘ni_n2-] 
X (1 - n3-i-)(l - n4-)vij{vij - V2j) 


^P4) 

(26) 


with magnetic field h = (pp — p-)/2 and normalization 
constant 


Cj_ — J ’^p-)- ( 27 ) 

The T matrix T{pi,P 2 ) describes s-wave scattering be¬ 
tween particle (pi,+) and particle (p 2 ,—)• In order to 
derive explicit expressions, the T matrix for ultracold 
fermions with s-wave contact interactions has to be used. 
The two-body T matrix 


To{pi,P2) 


AiThP 1 
TO a~^ +ik 


(28) 


is given in terms of the relative wavenumber k = \pi — 
P 2 \/ 2 h and the s-wave scattering length a. At weak cou¬ 
pling \a\ -)> 0 , To{pi,P 2 ) 'iTrh'^a/ TO is the regularized 
bare contact interaction. The BCS-BEC crossover goes 
from the weakly attractive BCS regime {l/kpa < —1) via 


the unitary Fermi gas (UFG, l/kpa = 0) to the repul¬ 
sive fermion branch [l/kpa > 1) above the BEG ground 
state [l3 | . At strong coupling the many-body T matrix is 
needed, which includes medium scattering. In the ladder 
approximation, which is the leading order of the large-A^ 
expansion in the number of fermion flavors [13,113 , the 
full T matrix reads 


dPp 


^p,+ '^p+p\+p2,— 


(27rfi)3 £p_^ -\- Sp,^ Sp ep-|_pj-|_p2 -I- 




(29) 


The T matrix is computed numerically, and Fig. [T] 
shows the resulting spin diffusion rate Tp/pp) in units 
of the Fermi frequency l/rp = Ep/h as the solid (red) 
line: it is highest in the trap center and decreases pro¬ 
portional to the density in the outer regions. At unitar- 
ity, the scattering rate is about 10 times larger than at 
weak coupling l/kpa = —2. In the Boltzmann regime, 

k/TD = [l - Psb - (/3eB)^exp(/3eB)Ei(-/3eB)] 

where Ei(a:) is the exponential integral and £b = k?/ma^ 
0 . At unitarity the scattering cross section decreases 
with temperature, and h/rp, = 4 -\/ 2 nA 3 /( 37 r/ 3 ) [13]. 

Third, the mean field Fimf = Wm{x,t)/h describes 
the precession of the spin current around the local mag¬ 
netization and is given by a momentum average of the 
real part of the many-body T matrix T{pi,P 2 ) over the 
momentum states between the majority and the minority 
Fermi surfaces, weighted by the velocity squared 0 , 


W = 


1 f d^pi d?p2 , . 

CTH J ( 2 ,(i)^ ("I, 

X (n 2 +- n 2 _)ReT(pi,P 2 )- (30) 


At weak coupling W = To(0,0) = dpti^a/m agrees with 
the bare interaction, which is real. At unitarity 1/a —>■ 0, 
To becomes purely imaginary and W appears to vanish 
along with ReTo. This is indeed observed at high tem¬ 
peratures, but at low temperatures the many-body T ma¬ 
trix T acquires a real part due to medium scattering, and 
kF ^ 0 0 . This is shown as the dash-dotted (blue) curve 
in Fig. m for weak coupling, W ~ ^pt/d'a/m is constant 
independent of the density and position in the trap, while 
at unitarity it decreases with density in the outer regions 
of the trap. Note that the very similar values for W at 
the trap center are coincidental: at weak coupling W is 
given essentially by the bare coupling, while at unitarity 
it is purely a many-body effect. 

Spin rotation is characterized by the dimensionless 
Leggett-Rice parameter 


ToWn 
1 = pn = --— 


(31) 


which is plotted as the dashed (magenta) line in Fig. [T] 
At weak coupling, 


1 


td oc 


2j.j4/3 ’ 


W 


a^n 


7 oc 


an 


1/3 ’ 


(32) 
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hence 7 becomes large and only weakly dependent on the 
density. At unitarity, on the other hand, 7 is a purely 
many-body effect, much smaller, and roughly propor¬ 
tional to the density. 


III. ANALYTICAL SOLUTIONS IN LIMITING 
CASES 


The evolution equations (m-dni) leave the density 
profile largely invariant, but the spin distribution changes 
dramatically as the magnetization is wound up into a he¬ 
lix according to the equations 

dtm + ^ Vj Jj -I- ax^z x m = 0 ( 33 ) 

3 

dtJj + a||P||Vj-m -I- a±P±Vjm + ujjXjm 

+ {axsz + —m) X Jj = -( 34 ) 

« Td 

The full numerical solution of these equations is pre¬ 
sented in Sec. HYi In order to gain a qualitative under¬ 
standing of the spin evolution in a trapped gas, it is 
instructive to consider first the approximate analytical 
solutions in the homogeneous case. 

If the scattering time T]j is much shorter than any 
other relevant time scale, for instance, the dephasing time 
T = (ZI-*-a^)“^/^ due to the helix, the current reaches a 
steady state and its time derivative vanishes in the ro¬ 
tating frame, dtJj = —ax^z x Jj. Defining = a±To 
and /im = —TuPtmf = iM one finds 

Jj -I- DQVjVn TD^^jXjm — fim x Jj = 0 . ( 35 ) 

This equation is solved by the steady-state current 


J, =- 




jm + fj.m x 


1 -I- jPm"^ 

-I-/rm(/im • Vjm)| — (36) 


where , /r, and td may still depend on position; the 
last term arises due to the trapping potential. Inserting 
this current into the continuity equation for the magne¬ 
tization (ESI) yields 


with effective diffusivity = Dq /{I + fj^m^). Since 
is constant in space, this is now a linear diffusion 
equation, albeit with a complex diffusion coefficient. 

In the spin-echo protocol, the gradient a winds up 
the transverse magnetization m_|_ = rrix + iruy into a 
helix along the X3 direction. A transverse spin current 
J3+ ^ 83171+ appears, which tends to smooth the helix. 
If J3+ has a component perpendicular to the local mag¬ 
netization m Ri rrizZ (at small tipping angle) it precesses 
around it at frequency = WrUz/h. At time a tt 
pulse around the y axis in Bloch space is applied; this is 
equivalent to reversing the sign of a. The subsequent 
time evolution unwinds the helix until the echo time 
te = where the transverse magnetization is again 
homogeneous. In the absence of spin rotation, 7 = 0, the 
modulus of the transverse magnetization A(t) = \m+{t)\ 
at time te decays as a cubic exponential @ 

A{te) = Aq exp(^- 

where Aq = A( 0 ). This result is approximately correct 
even for 7 7^ 0 if the tipping angle 0 is small, |m+| <C 
Imzl, and for short times. For finite 7, the magnetization 
decay is slowed down and is given by 0 


A{te) = AqJ-W ^77 exp 




6(l-f 72M2)_ 


Hte) = iMz In 


A{te) 

^0 


(39) 

( 40 ) 


where 77 = 72(Ao/7i)2/(l -\- y^M^) for polarization Mz, 
and yV(z) is the Lambert-IF function. This solution for 
the homogeneous system is used in the analysis of the ex¬ 
perimental data to fit the diffusivity Dq and the apparent 
Leggett-Rice parameter 7 from the measured magnetiza¬ 
tion decay and phase shift (j) of the trapped system. In the 
following the full spin evolution in the trap is computed 
explicitly to determine to what extent the homogeneous 
solution is still applicable to the trapped gas. 


IV. RESULTS 


dtm = —ax3Z x m + ’y^uj’^VjiTpXjm) 


Di 


J 

jm+yimxV jm+yLm{yLm-V jm)^. 


+V ,■ „ „ 

1 -I- y?m^ 

The Leggett solution @ is recovered in the homogeneous 
limit ijjj = 0, where remains constant in space: 


dtm = —axpz x m 




1 -I- ji^m- 


■{V^m -I- y,m x V^m}. 


In this case the longitudinal magnetization ruz remains 
unchanged, while the transverse magnetization m+ = 
nix + imy evolves as 


8 tm+ = —iax 3 m+ + L^^(l + ii 3 mz)'V^m+ ( 37 ) 


In the trapped gas, the assumptions which led to the 
analytical solutions in Sec. UTI] (steady-state current, ho¬ 
mogeneity) are not satisfied. Instead, a full numerical 
solution of the spin evolution equations, dMl) and (IMl) . 
is necessa^. The initial condition in the experimental 
protocol is a fully polarized cloud of fermionic atoms 
with a thermal profile, (HU) and ([ 22 ]), for tipping angle 9 . 
This distribution is stationary in the absence of an exter¬ 
nal magnetic field gradient a. Due to the density profile, 
m, also the coefficients aj_(a:), t±{x), and W{x) de¬ 
pend on the position in the trap. For a small tipping 
angle \m+\ <C \niz\ the gas remains almost fully polar¬ 
ized, hence the density profile and the coefficients are 
time independent even in the presence of the gradient a. 
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FIG. 2. (Color online) Transverse magnetization profiles on 
the BCS side Ijkpa = —2 at temperature {T/TF)i = 0.2 
and times tir just before the tt pulse (left column), and at 
the echo time ts = 2t-n (right column). Upper panels display 
the amplitude A = |m+|; lower panels, the phase angle cji = 
argm+. Solid curves in each panel are for a trapped gas at 
te = 0.2 ms to te = 1.2 ms, in comparison with a homogeneous 
gas at te = 1 . 2 ms (dashed line). 


In the numerical solution, the experimental parameters 
0 are used: gradient a = 1.67/rm“^kHz, trap frequency 
UJ 3 = 2tt X 750 Hz, a; = 27r x 470 Hz, roughly N ^ 40000 
atoms of with a Fermi energy of Ep = 2Trti x 29 kHz, 
and Thomas-Fermi radius Rtf = 5.1/rm. Even though 
the Fermi gas is initially fully polarized, the Fermi wave 
vector kp = 12fim~^ and the Fermi time Tp = 0.0087ms 
are taken for a hypothetical balanced gas of the same 
atom number in order to recover the standard relation 
n = kp/3n‘^ for the total density 0- At the lowest ex¬ 
perimental temperature, T Ri 280 nK, the initial reduced 
temperature is (T/Tp)i = 0.2 and A = 0.52/xm. The 
phase-space density at the trap center reaches A^n = 7.3, 
well within the quantum degenerate regime. 

In Fig. [21 the helical state of the local magnetization 
in the trap at 1/kpa = —2 is shown for different times 
ijr (left column). At this time, the tt pulse is applied 
(equivalently, the sign of a is reversed), and the subse¬ 
quent time evolution unwinds the helix until the echo 
time te = 2 < 7 r (right column), which is chosen between 
te — 0.2 ms and te = 1.2 ms. For comparison, the den¬ 
sity profiles for a homogeneous gas with te = 1.2 ms 
are shown as dashed lines. At time the (j) profiles 
of the trapped and the homogeneous gas have a simi¬ 
lar slope and are shifted toward negative phase angles 
by spin rotation. At time te, instead, the phase angle 
of the homogeneous gas is again homogeneous, but m+ 
for the trapped gas remains in a helical state. There¬ 
fore, the trap-averaged transverse magnetization decays 
quickly even though the local m+ is still sizable, and the 
apparent slope (j){te)/Mz\n{A{te)/AQ) is lower than in 
the homogeneous case, where it reaches the microscopic 







te [ms] 


ln(A/Ao)/7i 


FIG. 3. (Golor online) Amplitude and phase of the trans¬ 
verse magnetization m+ at temperature {T/TF)i = 0.2 and 
small tipping angle 6 = tt/SO (M^ = 0.995). Left: analyt¬ 
ical solution, (1391) (solid line), with Dq adjusted to fit the 
trap calculation (circles). Right: fit of 7 from the slope of 
the phase shift tp in Eq. (HDl) . Top, Ijkpa = —2 (BGS side); 
center, Ijkpa = 0 (unitarity); bottom, Ijkpa = 2 (repulsive 
branch). 


value 7 = —Wra^Tpijh according to Eq. (i40l) . As we see 
below, this effect leads to a saturation of the apparent 7 
for weak coupling. 

Note that the magnetization profiles are not symmetric 
with respect to X 3 \ this is due to the spin-rotation term 
in Eq. (1571) producing an additional phase shift on top 
of the gradient a. Consequently, if one reverses the sign 
of the scattering length a, the sign of 7 is approximately 
reversed, and the resulting magnetization profile is the 
mirror image with 2:3 1 —>■ —X 3 . Note also that at strong 
coupling \jkpa = 0 , the phase profile at te is nearly 
homogeneous at the trap center (not shown), and the 
kinetic coefficients for the trapped gas approach those in 
the homogeneous case. 

Figure |5| shows the decay of the trap-averaged trans¬ 
verse magnetization with time (left column), and the 
growth of the phase angle (p with the slope of the magneti¬ 
zation decay (right column). From these plots, the appar¬ 
ent Leggett-Rice parameter 7 = (j){te)jMz\n{A{te)jAo) 





















7 




to the homogeneous case (see below). Figure|4]also shows 
that the Leggett-Rice parameter changes sign, reflecting 
the sign of the effective interaction between quasiparti¬ 
cles which is attractive on the BCS side and repulsive 
on the repulsive branch [^. The sign change occurs for 
slightly positive values of I/Zcfo. 

The center panel in Fig. S] shows the bare transverse 
diffusivity . Again, the trap diffusivity is significantly 
lower than the homogeneous value (solid black line) at 
weak coupling and agrees with experiment. Once 7 is 
known, the bare diffusivity is found from the fit of 
the analytical solution (15^ to the magnetization decay 
in Fig. [3 

The bottom panel in Fig. |3 displays the ratio 

^ ?i 7 Wn 

—-r — x-■ (41) 

2mDQ 2ma± 



1 /kpS 

FIG. 4. (Color online) Interaction effect on spin transport: 
Leggett-Rice parameter 7 , bare diffusivity Dq, and interac¬ 
tion parameter Ao as a function of Ijkpa at degenerate tem¬ 
perature {T/TF)i = 0.2. Results for small {d = tt/SO; circles) 
and large {0 = rr/S; triangles) tippi^ angles are shown and 
compared to the experimental data [j (squares). 


can be read off as the slope of the curve through the 
origin. This value of 7 is then used to fit the magneti¬ 
zation decay on the left side of Fig. [3] to the analytical 
decay function (l39l) with diffusivity Dq. Here, I follow 
the analysis of the experimental data in Ref. 0 and fit the 
decay of the trapped gas to the homogeneous solution. 

In Fig. m the apparent Leggett-Rice parameter 7 ex¬ 
tracted above is plotted versus the interaction strength, 
from the BCS regime (left) via unitarity {l/kpa = 0) to 
the repulsive branch (right). For a trapped gas, I find 
that the apparent 7 saturates for weak coupling on both 
the BCS and the repulsive sides. This behavior differs 
qualitatively from the homogeneous case, where 7 con¬ 
tinues to increase linearly at weak coupling [cf. Eq. (1331) ] 
(solid black line). The trap calculation thus explains the 
saturation of 7 observed in experiments [J . The absolute 
value of 7 for the trapped gas depends on the tipping an¬ 
gle 9 of the initial polarization: for a large tipping angle 
0 = 7 r /3 as in the experiment, the Leggett-Rice parame¬ 
ter 7 saturates at larger values than in the case of a small 
tipping angle 9 = 7r/30. The value to which 7 saturates 
also depends on the strength of the gradient a, in contrast 


It measures the strength of the effective interaction irre¬ 
spective of the scattering time td and follows the sign 
change of 7 since Dq > 0. Again, the trap calculation 
agrees with the experimental data except in the instabil¬ 
ity region 0 < Ijkpa < 1.3 Q, while the homogeneous 
gas has more pronounced interaction effects. 

Spin rotation in a trapped gas is determined qual¬ 
itatively by the ratio between different length scales: 
the trap size Rtf ~ 5^m, the helix pitch ^heiix = 
211 j {ate/2) K, 4/im for te = 2 ms, and the mean free 
path £nifp = vpTD, which ranges from 0.5/rm at unitarity 
to 3.8/rm at l/kpci = —2. In simulations with a gradi¬ 
ent a larger than the experimental value, several helix 
pitches fit into the trap, Rtf ^ ^heiix, and the trap av¬ 
eraged diffusivity and 7 are equal to their homogeneous 
values as long as also £mfp < ^heiixj but saturate for larger 
^mfp- This is in marked contrast to longitudinal spin dif¬ 
fusion, where a scaling factor of about 5 was found to 
relate the trap-averaged and homogeneous diffusivity Ry 
[TtL m, [ 25 . For weaker gradients where t'heiix > Rtf, 
less than one helix pitch fits into the trap and the homo¬ 
geneous solution is reached not even at unitarity where 
Rtf ^ ^mfp- 

Finally, Fig. [3 shows the temperature dependence of 
the Leggett-Rice parameter at unitarity. 7 decreases with 
temperature, which is understood as follows: the two- 
body scattering amplitude is purely imaginary at uni¬ 
tarity and would imply a vanishing 7 ^ ReT. Hence, 
the observed finite value of 7 in the quantum degenerate 
regime i/TlTp)i < I is a many-body effect due to medium 
scattering. The presence of the medium enhances both 
the dissipative and the reactive effects of scattering at low 
temperatures, and more so in the homogeneous gas than 
in the trap, where only the core is strongly interacting. 

Similarly, the diffusivity reaches values of Dq ~ 2h/m 
in the quantum degenerate regime (center panel), in 
agreement with experiment, while the ratio Aq becomes 
large only at the lowest temperatures (lower panel). 






















(T/TF)i 


FIG. 5. (Color online) Temperature effect on spin transport at 
unitarity: Leggett-Rice parameter 7 , bare diffusivity Dq , and 
interaction parameter Ao as a function of the initial temper¬ 
ature (T/TF)i- At unitarity the results for small (9 = tt/SO; 
circles) and large {9 = tt/S; triangles) tipping angles differ 
very little; (blue) squares are experimental data [J. 

V. CONCLUSION 

To summarize, a spin-echo sequence in a trapped Fermi 
gas is modeled by a nonlinear and complex diffusion equa¬ 
tion for the transverse magnetization. The spin evolution 
exhibits the Leggett-Rice effect of strength 7 = /in, which 
appears to saturate for weak coupling, and a bare diffu¬ 


sivity much lower than expected for a homogeneous gas 
of the same temperature and interaction strength. These 
results are obtained without any fit parameters and agree 
very well with the weak-coupling data measured recen^ 
for a trapped gas of ultracold fermionic atoms Q. 
The present calculation provides an intuitive interpre¬ 
tation of the observed saturation of 7 : while the spin 
helix in the homogeneous gas is completely unwound at 
the echo time te-, the trapped gas remains partially in a 
helical state, with the average transverse magnetization 
A{te) strongly reduced and a smaller phase shift 
At weak coupling l/|A:i;’a| > 2 the kinetic theory em¬ 
ployed in this study is well controlled and includes the 
relevant interaction effects. 

At strong coupling l/kpa = 0 (unitarity), the Leggett- 
Rice effect is absent at the two-body level and arises only 
due to many-body medium scattering. The trap calcula¬ 
tion at strong coupling agrees with the experimental re¬ 
sults qualitatively; the remaining differences mw be due 
to reheating in a demagnetized Fermi gas d, m or due 
to interaction corrections to the equation of state which 
were not included in this study; they remain a topic for 
future work. 

Ultracold Fermi gases are thus ideal systems to study 
the interaction and temperature dependence of the 
Leggett-Rice effect and the spin transport coefficients. 
The spin dynamics of the trapped gas (Fig. [5]) is much 
more complex than in the homogeneous case and requires 
a numerical solution. The spin evolution simplifies near 
unitarity if Rtf > ^heiix > ^mfp, in which case the ho¬ 
mogeneous solution is recovered without any trap-related 
scaling factors. The decay of the trap-averaged magne¬ 
tization can be analyzed by the methods developed for 
the homogeneous case: the magnetization decay fits the 
homogeneous solution surprisingly well (Fig. [J]). This 
study shows that the extracted transport coefficients for 
the trapped gas can differ markedly, especially at weak 
coupling, from those of the corresponding homogeneous 
system (Fig. 0]). Conveniently for the interpretation, the 
sign of 7 , which reveals the repulsive or attractive char¬ 
acter of the effective interaction, does not change with 
the trap average. 

The author wishes to thank E. Taylor, J. Thywissen, 
S. Trotzky, and S. Zhang for stimulating discussions. 
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